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ABSTRACT 


The behavior of thin square perforated plates under 
the action of uniform shear deformation is studied experi- 
mentally and analytically using finite element analysis. 
Elastic Shear buckling strength is established as a function 
of the diameter of a round, centrally located hole in the 
plate. Post buckling behavior and the behavior of perfor- 
ated plates with various ring stiffeners are also studied 
experimentally . 
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1 . INTRODUCTION 

Thin plates occur In a wide variety of structures. 

While a good deal of our early knowledge about the behavior 
of thin plates has evolved from research directed to civil 
engineering structures such as plate girders, new applica- 
tions in aerospace structures and pressure vessels have ac- 
celerated research in recent years. 

Whatever the application, plate elements often block 
access to some other portion of the parent structure. Per- 
forations are made in plates to allow permanent passage of 
electrical, hydraulic, or other conduits, or simply to per- 
mit easy access for painting or servicing. Where a reduc- 
tion in strength is permissible, a hole might be placed in 
a plate to reduce weight. This is common in aerospace 
structures. In some cases, the reduction in the strength 
of the plate can be offset by the presence of a stiffener. 

Plate elements are often subjected to complex loading 
conditions. To understand the behavior of perforated plates 
under complex loads, we must first know something of their 
behavior under more simple loading conditions. The behavior 
of thin, square, perforated plates subjected to in-plane 
shear is studied experimentally and analytically in this 
report. The plates have clamped edges and centrally located, 
circular perforations. An experimental study of stiffened 
perforated plates is also presented here. 

A key aspect in the behavior of any thin plate is its 
buckling strength. Thin plates commonly become unstable 
at stress levels far below the proportional limit of the 
material. In some cases, the large transve -se deformations 
associated with buckling can be tolerated an I only the ulti- 
mate strength of plate is of any practical consequence. 
Elsewhere, aerodynamic or aesthetic considerations dictate 
that buckling cannot be allowed. 
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Plate buckling strength is usually expressed in terms 

of the critical stress, o or t . The critical stress is 

cr cr 

the stress which, when applied uniformly along the edges of 
the plate, produces instability. For an isotropic square plate of 
thickness t, width b, elastic modulus E, and Poisson's ratio 
v, this stress can be expressed as 


cr* cr 


= k 


« 2 Et 2 
I2(l-v 2 )b 2 


( 1 ) 


where k is known as the stability coefficient. Stability 
coefficients for compression and shear buckling of plates 
with various aspect ratios and edge conditions have been 
calculated and can be found in Reference 1. For the case 
of shear buckling of a square plate with clamped edges, 
this coefficient has been determined by Budiansky and Connor 
(Ref. 2) to be 14.71. 

Perforated plates, however, present some analytical 
problems. The case of the square perforated plate under 
edge compression has been studied quite extensively since 
1947. The most important works in this period are those of 
Levy, et al. (Ref. 3), Schlac* ,Ref. 4 and 5), Kawai and 
Ohtsubo (Ref. 6), and Vann and Vos (Ref. 7). The analytical 
methods used in these papers range from energy methods and 
numerical integration in the report of Levy, et al. to 
finite element analysis as used by Vann and Vos. 

The above work excluded consideration of shear buckling. 
Infinitely long shear webs have received some attention, but 
their relevance to the square plate is minimal and they will 
not be discussed in this report. The first serious attempt 
at predicting the critical shear stress of square perforated 
plates was made by Kroll (Ref. 8) in 1949. Using Timoshenko's 
energy method, which calls for the evaluation of two double surface 
integrals, Kroll obtained theoretical buckling coefficients for 
simply supported square plates with circular holes having 
diameters (d) equal to 1/8 and 1/4 the length of a side. 
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Numerical integration was used, and 22 points were 
needed in each octant to reduce the error for the known 
unpierced case to less than 5$. Analysis of the perforated 
plates indicated that the smaller hole reduced the coef- 
ficient by only 0.2$, while the larger hole caused a 22.6$ 
reduction in the coefficient. Hole reinforcements were 
also analyzed and were found to increase the coefficient 
by as much as 33^ $ over that for an unreinforced hole. 

In 1967, Rockey, Anderson, and Cheung (Ref. 9) solved 
the same problem using finite element analysis. With an 
idealization containing 56 triangular elements in the ana- 
lyzed quadrant, the solution for the unperforated case dif- 
fered from classical theory by 9 - 7% - Solutions for the per- 
forated cases differed markedly from those obtained by Kroll. 

Pour hole sizes were analyzed with simply supported and 
clamped edges. For the simply supported plate with a hole 
diameter to plate width ratio d/b = 0.125, Rockey, et al. 
found a reduction in the buckling coefficient of about 15$, 
whereas Kroll reported only a 0.2$ reduction. The two papers 
were in somewhat better agreement with respect to perforated 
plates with a hole diameter to plate width ratio d/b = 0.250, 
for which Rockey, et al. observed a 28$ reduction in the buck- 
ling coefficient, compared to Kroll 1 s 22.6$. 

Rockey, et al, investigated hole reinforcements, like 
Kroll, and found them extremely effective in canceling losses 
in buckling strength due to perforations. An experimental 
program was conducted and agreement with their analysis was 
good for the smaller diameter holes. Agreement was not as 
good for the larger holes. 

Both Kroll and Rockey, et al, studied uniform applied 
shear stress rather than the case encountered in testing, 
that of uniform shear deformation. This was most graphical- 
ly apparent in the investigation of Rockey, et al., wherein 
uniform stress was analyzed and uniform deformation tested. 
While this distinction has little effect on the classical 
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case of the unperforated plate, it has definite consequences 
when a hole is introduced. As the diameter of the hole in- 
creases, the stiffness of the plate continues to change along 
the edge and the two cases become more dissimilar. 

This distinction was observed for the perforated plate 
in compression by Schlack in 1964. It was taken into ac- 
count in his analysis, which used the Rayleigh-Ritz energy 
method, and was borne out by his testing. 

The study 'reported here is restricted to perforated 
plates subjected to uniform shear deformation. The choice 
of test specimen geometries and the experimental program 
are described in Section 2. An analytical approach to 
elastic stability of unstiffened perforated plates is pre- 
sented in Section 3, along with a summary of computed re- 
sults. The experimental results for both unstiffened and 
stiffened plate? ?»ro given in Section 4 and comparisons 
between analysis and test results are made. The final sec- 
tion of the report contains the concluding remarks. 

2. TEST SPECIMENS AND EXPERIMENTAL PROGRAM 
Specimen Size and Materials: 

The experimental program was cone .cted on 10” square 
aluminum plates with a thickness of 0.063”. The 10” size 
was chosen for convenience. It was felt that plates of 
this siz;e would be large enough for careful observation and 
instrumentation, and yet not too large to test in a standard 
testing machine. Larger plates would also have necessitated 
a larger and more costly test fixture. 

The material and thickness of the plates were chosen 
to ensure elastic buckling. Since very thin plates would 
be difficult to handle and would probably have initial im- 
perfections of more significance than thicker plates, the 
latter were considered more desirable than the former. 

Since buckling stress is proportional to the square of the 
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plate thickness and to the elastic modulus, thicker plates 
could only be used with lower modulus materials. Thus, 
aluminum, with an elastic modulus approximately one third 
that of steel, was chosen. 

Using Eq. 1, it was found that aluminum plates with a 
thickness of 0.063” and a side length of 10” would have a 
critical stress of 5*4 ksi. This is well below the yield 
stress and proportional limit of most of the common struc- 
tural aluminum alloys. 

The specific alloy chosen was 6061-T6, which has a 
relatively high yield point, an extensive linear post-yield 
region in its stress-strain curve, and a relatively small 
difference between its yield stress and ultimate strength. 
Hence, it approaches the elastic-plastic properties of 
structural steel and was chosen on that basis. 

Six tensile coupons were tested to evaluate the signifi- 
cant mechanical properties of the material. The properties 
of interest are the initial modulus of elasticity E, Poisson's 
ratio v, the proportional limit o » the yield stress a , and 

r 

the ultimate strength o u - The proportional limit was taken 
as the point at which tensile strain first deviated from 
linearity by 0.0001 and the 0.2!? offset approach was used to 
define the yield point. The above properties, which did not 
vary significantly with the direction in which the material 
was rolled, were E * 10.2 x 10^ p.s.i., v = 0.3> Op * 34,000 
p.s.i., o = 30,000 p.s.i. and o 0 45,000 p.s.i. 

J ** 

Three groups of specimens were tested. The first group 
were simply unperforated plates, (Figure 1, d = 0), used to 
establish "benchmark” values for the buckling loads. The 
second group of specimens was used to study the influence of 
unstiffened hole size on shear buckling capacity. Central 
holes with diameters (d) of 2, 4, and 6 in. were cut in these 
specimens, giving perforated plates with d/b ratios of 0.2, 
0.4, and 0.6. 
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The third croup of specimens (Pig. 2) consisted of 
perforated plates with stiffening material added around the 
perforations for the purpose of increasing the buckling load 
up to or beyond the buckling load of the unperforated plates. 
Five stiffened plates with four-inch holes were tested. All 
five stiffeners were flat rings made of the same material as 
the basic plates. 


Choice of Stiffeners 

There are several rules encountered in practice which 
can be used to determine the size of the required stiffener. 
One such rule is to replace the cross sectional area of the 
hole with a stiffener of equal cross section. This notion 
would have a theoretical foundation if buckling were entirely 
an axial phenomenon, but since the flexural rigidity of the 
plate dictates the buckling load, this rule lacks proper 
theoretical justification. 

Two approaches to the equal cross section rule were 
tested. Stiffener 1 was a 2" wide circular ring fastened 
to one side of the plate only, while stiffener 2 consisted 
of 1" wide circular rings fastened on both sides. Stiffener 
1 used 20 % more material than stiffener 2. 

Another common criterion is to replace the volume of 
material removed by the hole with a stiffener of equal 
volume. This rule has no more theoretical justification 
than the first, but it does result in a plate made of the 
same amount of material as the unperforated plate. 

Three different equal volume stiffeners were tested. 
Since all stiffeners were fabricated from material of the 
same thickness as the test plates, equal volume was equiv- 
alent to equal sur r ac«- area. Stiffener 3 was a simple 
circular ring, similar to stiffener 1 but only 0.828” wide. 
Stiffeners and 5 were fabricated from 5.012” squares to 
form square outer and circular inner boundaries. 
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Stiffener 4 was oriented with its edges parallel to 
those of the test plate, while stiffener 5 was oriented at 
45 degrees. Since shear is not a radial phenomenon, there 
is no reason to believe that a circular ring is the most 
efficient use of available material. Moreover, square 
stiffeners are easier to fabricate than circular ones. 

The reinforcing rings were uttachod to the plates with 
epoxy (Epon 907, produced by the Miller-Stephenson Chemical 
Company). In order to keep the thickness of the epoxy as 
uniform as possible, thin steel fibers ( 0 . 002 " dia.) were 
placed between plate and stiffener at clamping points. 

In each specimen, the epoxy was strong enough to hold 
the ring to the plate beyond the buckling load. However, 
the epoxy always failed in a brittle fashion shortly there- 
after, making it impossible to evaluate ultimate loads. 

Loading Apparatus and Test Procedures: 

The fixture used for testing the plates was designed 
to simulate the action r> f uniform shear deformation for 
plates with clamped edges. Such ideal behavior required 
the somewhat comple ' fixture shown in Fig. 3 . Specifically, 
the external edge members had to form a mechanism with 
hinge, located at the four corners of the plate to allow 
the plate to distort from a square shape to a diamond shape 
as shear load was applied. 

The desired loading and edge restraint was accomplished 
by using a frame that consisted of two separate pairs of 
frames. An exploded view of the system is shown in Figure 
4, including a typical test specimen. The test specimen 
had a 3 in. wide strip along each edge to enable bolting of the 
specimen to the frame. The corners of the xesultit.g 1 6 in. 
square plate were cut out, as shown in Figure 5, to eliminate 
any extraneous stresses that might arise from distortion 
of the areas outside the basic 10 in. dimension 0 * the plate 
as shearing Ttortions were introduced in the specimen. 
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The inner pair of frames were intended to force the 
specimen edges to displace linearly. They were made of 
3 ” x 3 / 4 " steel, formed the 10" square boundary of the 
test plates and provided the clamped edges over a 3" strip 
around the perimeter of the plates. Attachment of the plates 
was made by a total of 52 1/2" diameter bolts arranged in 
two rows. The holes for these bolts had to be located with 
extreme accuracy to assure proper fit between the fixture 
and the specimens. Moreover, the bolts had to fit into the 
holes with a minimum of play so that the plate could bear on 
the bolts in the event that frictional forces were not suf- 
ficient to transfer load into the plate. It was expected 
that possible slippage at the frame-plate interface would 
not happen until after the plate had buckled. 

The outer pair of frames, made of 2-1/2” x 1/2” steel, 
were connected to the inner pair through the central row of 
bolts. It was the outer pair chat contained the corner 
pivots in the form of 1" diameter hardened steel pins. These 
pins also provided the diagonal load points. 

The test set-up was assembled vertically in a deforma- 
tion controlled testing machine capable of exerting a tension 
load up to 30,000 pounds. 

The fittings on either end of the fixture were each 
fabricated from three separate pieces of steel so that they 
could easily be moved and assembled. With just one side of 
the fixture (inner and outer frames) assembled, the test- 
ing machine could be adjusted until the 52 connecting holes 
in the fixture and specimen lined up. The machine was kept 
in this position while assembly of the fixture was comt ?ted. 
There were no forces on the plate at this point and no 
further computations were necessary to determine the point 
of zero load. 

The test specimens were instrumented with 1” gage 
length strain gages placed at points of local strain maxima 
to monitor in-plane behavior and check for local yielding. 
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In certain cases, strain gages were placed back-to-back on 
opposite faces of the test plates to aid in the detection 
of buckling. The points of maximum expected strain were 
determined from the computer analysis. 

Dial gages were used to measure out-of-plane deflections. 
These gages were placed at locations of maximum transverse 
deflection, which were determined from the eigenvector ob- 
tained from the computer program. Strain and deflection 
gage readings were taken at small increments of tensile load. 

At each increment, the load was permitted to stabilize be- 
fore the readings were taken. The results obtained and 
their interpretation in terms of buckling loads will be de- 
scribed in Section 4. 

3. ANALYSIS 

Formulation Basis 

The analyses for initial buckling of the unperforated- 
and perforated-unstiffened plates were performed with use of 
a finite element computer program whose theoretical basis is 
described in References 11 and 12. We briefly outline in 
the following the key features of this theory and those fac- 
tors which are relevant to the subject analyses. 

The finite element used in this work is a triangular 
thin shell element, portrayed in Figure 6. The flat form of 
the element of course applies in the present work. The 
relationships for this element are constructed in the stiff- 
ness matrix format on the basis of a generalized potential 
energy approach. In accordance with this approach, an ele- 
ment stiffness matrix is first formed using rather simple 
displacement functions. These violate conditions of inter- 
element continuity of displacement when the element is 
Joined to others of the same type. The continuity is "restored” 
by writing equations which stipulate the required continuity 
on each of the element boundaries. Such equations are con- 
straint equations which are incorporated into the global 
analysis by means of the Lagrange multiplier method. 
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In this work the u, v and w displacement fields of 
the triangular element are chosen in the form of complete 
cubic polynomials. The vertex node points each have 9 
degrees-of- freedom (the translational displacement and the 
first derivatives with respect to the normal coordinate, 
for each of the three translational displacement components ) , 
and a centroidal coordinate with the values of the three 
translational displacements as degrees-of- freedom. This 
formulation lacks interelement continuity of normal angular 
displacement across the element boundaries. This continuity 
is restored by writing a constraint equation which requires 
that the angular displacements at the midpoints of the sides 
of adjacent elements be equal. 

The formulation for linear elastic instability analysis 
involves the formation of a linear stiffness matrix [K] , and 
a "first -order" geometric stiffness matrix [N^]. The geomet- 
tric stiffness matrix is a function of the applied loads and 
is calculated on the basis of a reference value of such loads. 
The intensity of the applied loads to cause instability is 
X times the reference value so that the value of the geometric 
stiffness matrix at buckling equals A[N^]. The condition for 
instability is that the determinant of the matrix [K] + A [N^] 
equal zero, i.e. 

| [K] + A [Nj] | = 0 (2) 

The computer program developed to implement the present 
element formulation seeks a zero determinant iteratively. 
Analysis is first performed for A * 0 in order to establish 
a basis for computing [N^], using a reference value of the 
applied loads. A is then set equal to 1.0 and the deter- 
minant is evaluated. These two points on the A vs. deter- 
minant plot are sufficient to start a Lagrangian interpolation 
scheme for predicting a zero determinant. Thus linear inter- 
polation is used first, followed by parabolic interpolation 
when three points are available, and so on until convergence 
Is reached. 
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The curve of load versus determinant is convex, which 
is unusual in buckling phenomena. In a compressive buckling 
problem, the curve would be concave. The difference is that 
in the case of shear, the direction of the load is of no 
consequence, while if a compressive load were changed to 
tension, buckling would never take place. For shear buckling, 
the determinant iteration curve must be symmetric about the 
vertical axis and, hence, must be convex. 

Finite Element Idealization 


The inherent symmetry of the problem, illustrated 
in Figure 1. simplifies the analysis somewhat. With the axes 
oriented along the diagonals as shown, each quadrant is iden- 
tical and only one need be analyzed. Care must be taken in 
the imposition of the proper boundary conditions along these 

axes. Along the x axis, for example, symmetry dictates that v 
dw 

and -r- must both be zero. Since v is zero along the entire 
length of the axis, must also be zero. Finally, there can 
be no shear along the x axis. Since 


xy 


= il 

3x 


3u 

ay 


(3) 


and — has already been determined to equal zero, must 
also be zero. Four similar boundary conditions along the y 
axis can be determined. 


An idealization using three intermediate radial lines 
dividing the plate edge into four equal lengths was used to 
simplify the introduction of various hole sizes. This is 
shown in Figure 8 for the four analytical models of this study, 
the unperforated plate and plates with three different hole 
diameters. In this way, the diameter could be changed from 
0 to 0.6b by merely scaling coordinates and without changing 
the basic style of the grid. 

The analysis was limited to holes no larger than 0.6b, 
in consequence of limits on the geometric form of individual 
finite elements. As triangular finite elements become elon- 
g .rid, their solution accuracy tends to deteriorate. For a 
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simple constant stress element, this effect becomes notice- 
able when the aspect ratio nears three but with the more 
sophisticated element used in this analysis this effect should 
not produce significant errors until the aspect ratio reaches 
five or six. The highest aspect ratio used in the present 
analysis was less than five. If a larger diameter hole were 
introduced, this difficulty would have to be dealt with and 
the basic style of the analysis grid would have to be changed. 

Uniform shear deformation was applied as prescribed u 
and v displacements at each of the five edge nodes. This 
represents a departure from the case of uniform shear stress, 
wherein x and y forces would be applied to these five nodes. 
Some insight into the source of the prescribed displacements 
can be gained by referring to Figure 7, which shows that the 
edge "slides" along the axes, remaining straight, and without 
changing length. 

For small deformations, the displacement of each point 
has the same component along the edge. Hence, 1-1* is equal 
to 5-5' ^nd the u and v displacements of each point conform 
to the ratios listed below. 


U 1 

* 

- 1.0 

V 1 = 

0.0 

u 2 

s 

- 0.75 

V 2 = 

0.25 

u 3 

X 

- 0.50 

v 3 

0.50 

u 4 

s 

- 0.25 

v 4 - 

0.75 

u 5 

s 

0.0 

V 5 = 

O 

• 

iH 


T 1 e nature of the eigenvalue iteration used in the computer 
program insures that these ratios do not change. Simulation 
cf clamped supports along this edge was accomplished by sup- 
pression of w, and -g— at each of the five nodes. 
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Convergence Study for Unperforated Plate 

The first objective in the analysis was to determine the 
convergence properties of the solution. As in any finite ele- 
ment analysis problem, it is essential to establish that the 
results were converging to the correct rolution as the grid 
was refined, and that the chosen grid refinement yields a 
sufficiently accurate answer. 

The unperforated plate was used for this purpose. For 
an unperforated plate, uniform shear stress causes uniform 
shear deformation. Therefore, the two cases are identical, 
and the buckling load can be computed using Eq. 1. 

Finite element solutions were obtained using gridworks 

of 12, 20, and 36-elements, respectively. Hie 12 and 20 

element meshes are shown as Figures 9a and 9b, respectively. 

The 36 element mesh appeared as Figure 8a. The classical 

"exact" solution, from Ref. 1, is T = 4,800 lb. The cal- 

cr 

culated results are T = 5,226 lb., S^l 1 * lb., and 4,954 lb., 
from the gridworks at Figs. 9a, 9b, and 8a, respectively. 

The results, shown graphically in Figure 10, demonstrate that 
the numerical solutions are indeed converging to the classical 
solution. The 36 element mesh produces an answer which is 
only 3 - 2 % high. 

Perforated Plate Solutions 

With sufficient convergence established, the buckling 
loads for the perforated plates could now be determined. 

The primary goal of these analyses was to establish the 
relationship between the buckling load and the hole diameter. 
To discuss these results, we must examine the manner of 
definition of critical loading. 

In the case of applied uniform shear stress, it is per- 
fectly reasonable to think in terms of critical shear stresses 
and buckling coefficients. There is a distinct value of 
shear stress which, when applied uniformly around the per- 
imeter of the plate, causes instability. 
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When the plate is under the action of uniform shear de- 
formation, shown in Figure 7, this is no longer the case. 

The shear stresses imposed by this deformation can be 
quite nonuniform. Moreover, they can be accompanied by 
direct stresses of the same order of magnitude. The crit- 
ical stress concept is inappropriate under these conditions. 

The critical displacement of one of the corners, or 
the critical rotation of an edge could accurately define 
the point of buckling, but neither of these values is of 
much interest to the design engineer. He is primarily 
concerned with the load which can be applied before the 
plate buckles. For this reason, the shear buckling strength 
will henceforth be expressed in terms of the critical diag- 
onal tension force, T . This will be defined as the total 

cr 

force which, when applied across either diagonal of a square 
plate which is constrained to uniform shear deformation, 
causes elastic instability. 

The relationship between buckling strength and hole 
diameter is shown in Fig. 11. This curve, particularly in 
the region of small diameters, is approximately parallel 
to and in fairly close agreement with the findings of Rockey, 
et. al. (Ref. 9), for uniform shear stress which are also shown 
in Figure 11. Although the percentage discrepancy of the 
two solutions grows with d/b ratio, the absolute difference 
between them remains essentially constant. In assessing 
this error, it must be recognized that different types of 
elements and gridwork refinements are used in the two ap- 
proaches and the associated finite element discretization 
error is likely to be on the order of the indicated dis- 
crepancy. Differences in the imposed boundary conditions 
might also play a role. 

The properties of Fig. 11 can be given an intuitive 
interpretation. A small hole would not be expected to 
greatly reduce the flexural stiffness of the plate, but 
would rather increase the stresses by the factor (1 -d/b). 
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This is precisely what is observed. As the diameter in- 
creases, the reduction in plate stiffness becomes more 
significant, causing the carve to become concave. 

The buckling load does not give all the pertinent 
data about the behavior of perforated plates. Also of 
interest is the in-plane stiffness, which can be thought 
of as the diagonal tension force per diagonal extension. 
Figure 12 shows that the prebuckling stiffness, defined 
as the ratic of the diagonal force T and the corresponding 
corner point displacement u, does not change much when a 
small hole is introduced, but decreases rapidly with in- 
creasing hole diameter. Hence, there is an inherent trade- 
off involved in using plates with large perforations. Al- 
though the loss in buckling strength is less than propor- 
tional to the d/b ratio, the increased flexibility might 
make such a plate undesirable. 

Further insight into the behavior of perforated plates 
can be gained through the examination of the buckling modes. 
The normal displacement contours for the four cases are 
shown in Figure 13* The oval nature of the buckling contours 
is readily verified by intuition. The region of negative 
deformation is, however, a bit more interesting. The fact 
that the area of negative deformation decreases in size as 
the hole diameter increases offers a clue to the region's 
behavior. It is a consequence of diagonal compression alonr 
the horizontal axis, which also decreases with increasing 
hole size. The unperforated plate deforms into three half- 
waves along the horizontal diagonal. When a hole is intro- 
duced, the interior edges are permitted to approach each 
other on the horizontal diagonal, reducing the buckling 
effect along the axis. 

4. RESULTS AND ANALYSIS-TEST COMPARISONS 

The test results on the perforated, unstiffened plates 
are presented in this section and compared with analytical 
values. This material is followed with a summary of the 
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experimental results for the stiffened perforated plates; 
no analytical results are available for the latter. 

Perforated, Unstiffened Plates: 

Buckling is a gradual process due to initial imperfec- 
tions. These initial imperfections were more pronounced and 
the buckling loads more difficult to pinpoint in the plates 
with the larger holes. In all cases, buckling could be de- 
tected through the strain gage readings long before it could 
be discerned visually. Typical plots of load against trans- 
verse deflection are shown for the four plate configurations 
in Figure 14. These plots indicate the gradual nature of 
buckling in each case. 

The problem of buckling point definition has been en- 
countered before and various criteria have been established. 
For example, Rockey, et al. (9), defined th buckling load 
as the intersection of the projections of the straight por- 
tions of the load vs. lateral deflection plot. This method 
has little justification but does produce an answer which is 
near the center of the transition range. Figure 15a, a load 
vs. lateral deflection plot for the unperforated plate, shows 
that this criterion would indicate that the buckling load for 
this plate was 5,200 pounds. 

Another criterion defines the buckling load as the in- 
flection point on the load vs. lateral deflection curve. 

(See Fig. 15b). This point does have physical meaning. It 
is the point of minimum plate rigidity in that the curve has 
a steeper slope on either side of it. Buckling points de- 
fined in this way are on the upper end of the transition 
range. Figure 15b shows that the buckling load for the un- 
perforated plate is 5,^00 pounds when this criterion is used. 

Strain gage data can also be used to define the buckling 
load. When two strain gages are placed back-to-back on op- 
posite faces of a plate in a region where buckling produces 
curvature, the gages indicate equal strains until the curva- 
ture appears. The strain reversal criterion defines the 
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buckling load as the point at which the strains read from 
one of the gages change direction. Figure 15c, which repre- 
sents strain gage data from the unperforated plate, indicates 
that this method would define the buckling load as 5,100 
pounds. One difficulty met with the strain reversal criterion, 
is placement of strain gages. They should be placed at the 
point of initiation of buckling-induced curvature, but this 
point is not always known in advance. 

Hence, the buckling load for the unperforated plate could 
be called 5,400, 5,200, or 5,100 pounds, depending on which 
criterion was used to define it. Although the inflection 
point method does not give the best answer for this case in 
terms of agreement with theory, it does have physical signif- 
icance and always gives a well-defined point at the top of 
the transition region. For these reasons, the inflection 
point method is used to define the buckling point in this 
report . 

The experimental versus analytical buckling loads for 
the four unstiffened plates were 5,400 lb. vs. 4,954 lb. for 
d/b * 0, 4,100 lb. vs. 3,940 lb. for d/b - 0.2, 2,700 lb. 
vs. 3»Ol4 lb. for d/b - 0.4, and 2,200 lb. vs. 2,482 lb. for 
d/b * 0.6. Figure 16 shows this information graphically. In 
view of the difficulties encountered in defining experimental 
buckling loads, the results are in good agreement. The max- 
imum deviation is about 10$. 

The strains encountered prior to buckling do not agree 
quite as well with analysis. Table 1 compares strains from 
a selected gage in each test with the corresponding strain 
from the computer analysis. In each case, the comparison is 
made for a centrally located gage at about one half the? 
buckling load of the plate. The maximum deviation in this 
comparison is 30?. This larger deviation than in other solution 
parameters (e.g. the critical structure) may be due to the 
presence of combined bending-flexure and the finite size of 
the gage, among other experimental effects. Also, in finite 



18 


element analysis with the displacement method, the strains 
are expected to be less accurate than the displacements, since 
they are given by the derivatives of the displacements. 

Buckling was the phenomenon of primary interest through- 
out the testing program. However, post buckling behavior 
was also of interest, and each plate was tested to failure. 
Figure 17 shows, sequentially, the behavior of the plate 
with a six inch hole through the onset of buckling to ulti- 
mate load and failure. The vertical ridge which can be seen 
at the end of the sequence was typical for the perforated 
plates. 

Figure 18 shows each test plate after being tested to 
its ultimate strength. Each perforated plate exhibited yield- 
ing in the direction of loading on either side of the hole. 
Eventually, the load decreased and fracture occurred at the 
same location. The unperforated plate yielded in bearing 
around the connecting bolts. Since this was an edge failure, 
it is difficult to draw conclusions on the ultimate strength 
of the unperforated plate. 

Each perforated plate exhibited an ultimate strength 
equal to about four times its buckling capacity. The var- 
iation of ultimate strength with d/b ratio is shown in Figure 
16 . 

In studies dealing with perforated plates under uniform 
shear stress, the relationship has been given as a plot of 
the buckling coefficient, k, against the d/b ratio. In this 
way, the results have been non-dimensionalized . In the 
present study, stresses are not constant and the buckling 
coefficient is not meaningful for reasons previously dis- 
cussed. The following parameter has therefore been defined: 


N 


T (l-v 2 )b 
cr 



2 

This dimensionless parameter is equivalent to k/?u /12 for 
the case of the unperforated plate. 
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The analytical results are plotted in this form in 
Figure 19* Also shown are the experimental results and the 
classical solution for the unperforated plate. In this 
figure, the measured thickness of each plate and the observed 
mechanical properties of the material were used in calcula- 
ting each experimental point. 

Using Figure 19, it is possible to evaluate the relative 
efficiencies of perforated and unperforated plates. It must 
be remembered that the introduction of a hole into a plate 
reduces both the buckling strength and the amount of material 
used. Consider, for example, a 100” square steel shear panel 
which must carry a diagonal tension load of 100,000 pounds. 
For an unperforated plate (d/b=0), the required thickness 
would be 0.257" and 2,570 cubic inches of material. 

if the same shear panel were fabricated with a 50" 
diameter hole, d/b=0.5 and N ■ 9.4. Hence, the required 
thickness would be 0.318" and 2,560 cubic inches of material 
would be needed. This represents a slight saving over the 
unperforated plate. This is not a significant saving, al- 
though perforated plates would appear more efficient at 
smaller b/'t ratios. However, other factors, such as stiff- 
ness, must be considered in the comparison. The in-plane 
stiffness of the perforated plate in the example would be 
much lower than that of the unperforated plate. Figure 12 
indicates that the introduction of such a perforation would 
lower the stiffness by about 40JK if it were not accompanied 
by a corresponding increase in thickness. Since the thick- 
ness increased by only 24£ when the hole was introduced, the 
redaction in stiffness would be quite significant. 

A similar reduction in the ultimate strength of the 
plate would be expected. Figure 16 indicates that the reduc- 
tion in ultimate strength would be even greater than the 
reduction in stiffness, but exact evaluation of this reduc- 
tion is not possible because of the edge failure encountered 
when testing the unperforated plate. 
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Moreover, initial imperfections are likely to be more 
serious in perforated plates. While diagonal tension tends 
to eliminate imperfections along the loaded diagonal of un- 
perforated plates, this effect is not present in perforated 
plates , and perforated plates would also be more likely to 
have large imperfections in the first place. Although initial 
imperfections do not change the buckling load, they have 
severe effects on the plate behavior prior to buckling, caus- 
ing large out-of-plane deformations. 

Stiffened Perforated Plates 

We consider now the experimental results for t tif- 
fened perforated plates, (Figure 2), for which no analysis 
results are available. Figure 20 gives the load vs. lateral 
deflection ( 6 ) of the unstiffened plate to those of the 
same plate reinforced with stiffeners 1 and 2. Lateral de- 
flections were measured at the expected point of maximum 
displacement. From this figure, it can be seen that both 
reinforcing schemes were quite effective in increasing the 
buckling loads of the perforated plates. 

Evaluation of the buckling load for stiffener 2 was 
difficult., but no more so than it had been for the unstif- 
fened plates. Due to the eccentricity produced by attaching 
t stiffener to only one side of the plate, however, the plate 
using stiffener 1 buckled quite gradually and determination of 
the inflection point involved a great deal of judgment. It 
must be noted that evaluation of the buckling loads for this 
and all other eccentrically stiffened plates carry with them 
an inherent Judgment or uncertainty factor on the order of 
15 %. 

Despite these difficulties, the inflection points were 
located and the buckling loads evaluated. It was found that 
stiffeners 1 and 2 increased the buckling loads by 237 ? and 
318?, respectively, over the same perforated plate with no 
stiffener. 



Figure 21 shows the load vs. lateral deflection curves 
for stiffeners 3> and 5, along with a similar plot for 

the unstiffened plate. As this figure indicates, the snape 
and orientation of the stiffener made little difference in 
behavior. Each stiffened plate exhibited the same gradual 
buckling that was seen in the plate using stiffener 1. 

Each increased the buckling load by about 60? so that each 
stiffened plate had a buckling load approximately equal to 
that of an unperforated plate. Table 2 lists the inflection 
point buckling loads for all five stiffened plates as well as 
the buckling strengths of unstiffened plates with and without 
the diameter hole. 

Although the numbers found in this table indicate that 
the circular ring is the most efficient equal volume stif- 
fener, the differences are too small to draw final conclu- 
sions. Indeed, if either strain reversal or the method of 
Ref. 9 had been used to define the buckling point, stiffener 
5 would have appeared to be the most efficient, followed by 
stiffeners 3 and respectively. 

5. CONCLUDING REMARKS 

The behavior of thin, square, perforated plates under 
the action of uniform shear deformation has been studied 
analytically and experimentally. The relationship between 
the elastic buckling strength and the diameter of a round, 
centrally located hole has been established using finite ele- 
ment analysis. Experimental results were in reasonable agree- 
ment with the analysis results. 

The experimental program on stiffened plates indicated 
that a ring of the same material and thickness as the plate 
can be extremely effective in canceling the loss of buckling 
strength which accompanies the introduction of a perforation. 
Rings made of the same amount of material as the removed hole 
can i *ase the buckling strength back to that of an unper- 
fora • ,;late and, in this study, alternate reinforcing 
schemes increased the buckling load by as much as 318?. 
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Whereas buckling of unstiffened plates involves only one 
variable, stiffener design Includes many, such as stif- 
fener size, shape, location, or material type. 

There is a considerable amount of additional research 
needed in the field of shear buckling of perforated plates. 
Problems that need attention include: 

a. analytical treatment of stiffened plates, with 
particular emphasis on optimum stiffener geometry. 

b. analytical and experimental study of plates with 
holes 01 differing shapes (such as rectangles) and at loca- 
tions other than the center of the plate. 

c. treatment of inelastic buckling and ultimate strength 
of plates. 

Experimental studies of shear buckling are quite expen- 
sive and time consuming. They must be supplemented as much 
as possible with analytical approaches. It appears that 
finite element analysis capabilities could play a major role 
in any future i ^search in plate buckling. 
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FIGURE 1 - Geometry of Test Specimens. 
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FIGURE 2 * Reinforcement Patterns Around Perforations in 
Plates with d/b ■ 0.4. 















FIGURE 3 - Testing Apparatus 
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FIGURE 6 - Plate Element Nodes and Degrees-of-Freedom 
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FIGURE 7 - Mode of Deformation of Specimen Due to 
Imposed Edge Displacements. 



a. Unperforated b d/b =0.2 



c. d/b = 0.4 d. d/b = 0.6 


FIGURr. 8 - Finite Element Gridworks 
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FIGURE 9 - Finite Element Gridworks for Convergence Study 
of Shear Buckling Calculations for Unperforated 
Plate. (See Fig. 8a for 36-Element Gridwork) . 
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FIGURE 10 - Comparison of Various Finite Element Models 
with Classical Solution. 
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FIGURE 19 - Comparison of Experimental and 
Analytical Buckling Loads for 
Various d/b Ratio. 
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KTOUHK 21 - T vs. 6 


for Equal Volume Stiffeners 


STRAIN 

LOCATION 



T, lb, 


2,500 376x10 


STRAIN 


350x10 


2,000 

278x10*^ 

1,500 

354x10*6 

1,200 

! 

473x10*6 


319x10 

462x10 

531x10 


TABLE 1 - Comparison of Analysis and Test Strains at 
Selected Gage Locations. 


STIFFENER 

NONE 

NONE 

#1 

#2 

#3 

#4 

#5 


T , lb. 
(EXPER.) 

5,400 

2,700 

9,100 

11,300 

5,000 

4,800 

4,500 


TABLE 7 - Inflection Point Buckling Loads for Both Unperforated 
Plate and Stiffened and Unstiffened Plates with 4 -Inch 
Perforation. 












